Series Expansion Analysis of a Frustrated Four-Spin- Tube 
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We study the magnetism of a frustrated four- leg spin- 1/2 ladder with transverse periodic boundary 
conditions: the frustrated four-spin tube (FFST). Using a combination of series expansion (SE), 
based on the continuous unitary transformation method and density-matrix renormalization group 
(DMRG) we analyze the ground-state, the one-, and the two-particle excitations in the regime of 
strong rung-coupling. We find several marked differences of the FFST with respect to standard 
two-leg ladders. First we show that frustration destabilizes the spin-gap phase of the FFST which 
is adiabatically connected to the limit of decoupled rung singlets, leading to a first order quantum 
phase transition at finite inter-rung coupling. Second, we show that apart from the well-know 
triplon branch of spin-ladders, the FFST sustains additional elementary excitations, including a 
singlon, and additional triplons. Finally we find, that in the two-particle sector the FFST exhibits 
collective (anti)bound states similar to two-leg ladders, however with a different ordering of the 
spin-quantum numbers. We show that frustration has significant impact on the FFST leading to a 
flattening of the ground-state energy landscape, a mass-enhancement of the excitations, and to a 
relative enhancement of the (anti)binding strength. Where possible we use DMRG to benchmark 
the findings from our SE calculations, showing excellent agreement. 
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PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Dg, 75.10.Kt, 02.70.Wz 



I. INTRODUCTION 

Ever since it has been realized that there are 'sur- 
prises' on the way from one- to two-dimensional quan- 
tum magnets^, spin ladders have attracted an enormous 
interest. Two-leg ladders, such as Sri4_ x Ca a; Cu2404i— , 
have been under intense scrutiny mainly because of the 
potential interplay between their spin-gaped ground state 
and the occurrence of superconductivity^. N-leg ladders 
with N > 2 are of particular relevance, not only because 
they allow for a generalized test of Haldane's conjecture^, 
but also because new, tube-like lattice structures can be 
realized if periodic transverse boundary conditions ap- 
ply, such as in [(CuCl 2 tachH) 3 Cl]Cl2^ and CsCrF^ with 
N = 3, Cu 2 Cl 4 -D 8 C 4 SO^ with N 



wnii iv = 4, and Na2V307«ii 
with N = 9. For N = 3, magnetic frustration surfaces 
as an additional ingredient, already for tubes with only 
nearest neighbor exchange, leading to a rich variety of 
phenomena not present in two-leg ladders^ - — . 

Experimentally, the four-spin tube Cu2Cl4-DgC4S02 
has been suggested to display frustrating antiferromag- 
netic next-nearest neighbor exchange^. Theoretically, 
unfrustrated four-spin tubes have been considered in two 
seminal paper a 15 : 16 , both, in the weak and the strong 
rung-coupling limit. Magnetic frustration, however, has 
not been considered in these studies. Therefore, in this 
work, we perform a first analysis of a frustrated four-spin 
tube (FFST) 
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with a lattice structure and exchange couplings ji m as 
shown in Fig|TJ Spin- 1/2 moments are located on the 
solid circles and all couplings, jo, j\ and ji are antiferro- 
magnetic (we set jo = 1 hereafter) . This FFST is simpler 
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Figure 1: Frustrated four-spin tube. Solid circles represent 
spin-1/2 moments. Plaquettes (bold black lines) are coupled 
by nearest (ji) and next nearest (j'2) antiferromagnetic ex- 
change, blue and red lines. On-plaquette couplings (jo) are 
set to unity. 



than the one proposed for Cu2Cl4-DgC4S0 2 , where only 
part of the surface squares experience diagonal exchange, 
and the leg-couplings seem to be in-equivalent^. Apart 
its relation to existing materials, the FFST is of interest 
as a ID variant of the anisotropic triangular lattice on 
the torus with four site circumference, i.e. cutting the 
tube longitudinally one obtains an anisotropic triangular 
lattice strip. 

For j'1.2 -C 1, the FFST resembles a chain of weakly 
coupled four-spin plaquettes each of which displays a 
singly degenerate singlet ground state, separated by a gap 
of jo = 1 from the first excited triplet. Therefore pertur- 
bation theory in j\^ applies. Motivated by this we inves- 
tigate the FFST by series expansion (SE) in j'1.2. More- 
over, we corroborate our approach and gauge our SE re- 
sults by employing density-matrix renormalization group 
(DMRG) calculations. The structure of the paper is as 
follows. In section [Hi we clarify the region of applicabil- 
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Figure 2: Phase diagram of the FFST. Solid lines: transitions 
at the classical level (all are discontinuous). Wave vectors 
A, B, C and D label pitch angle of classical phase. Circles: 
first-order quantum critical line from DMRG (see text). Inset: 
ground state energy versus ji from DMRG for the quantum 
model at ]2 = ji, showing first order transition at ji = jf « 1. 

ity of the SE. Section Mil details our SE method. Results 
are presented in section IIV| including the ground state 
energy of the FFST in sub-section IIV A[ the one-particle 
excitations in sub-section IIVB1 as well as two-particle 
states in IIV CI Conclusions are presented in section [Vj 
For completeness, technicalities of the two-particle SE 
calculations are included in the appendix IVHI 



II. COUPLED PLAQUETTE REGIME 

Proper application of SE hinges on the adiabatic renor- 
malization of the bare starting state in terms of the 
coupling constants. In case of intervening second order 
quantum phase transitions, the SE can be used directly 
to limit its range of applicability in terms of diverging 
susceptibilities or vanishing elementary excitation gaps. 
In case of a discontinuous or first order transitions, SE 
based on a single bare ground state fails to signal any 
transition. To put our SE on firm grounds a-priori, we 
therefore search for potential first-order quantum phase- 
transitions of the FFST at small ji 2. 

To this end, it is instructive to first consider the classi- 
cal phase diagram of the FFST. We allow the spin struc- 
ture to be a spiral, which, due to SU(2) symmetry can be 
considered to be planar S = 5(cos(Q-r), sin(Q-r), 0) with 
r = l x R x + lyRy, where R x>y = (1, 0), (0, 1), l x £ N, and 
l y = [1, . . .4]. The transverse pitch vector Q y has to be 
discretized according to (0, 1, 2, 3)ir/2. The ground state 
energy is e c g = cos(Q y ) + j ± cos(Q x ) + j 2 cos(Q x + Q y ). 
From this, four classical phases result, shown in Fig. 



[5] a (tt, 7t) antiferromagnet for 32 ^ (1 + 2ji)/(2(ji + 
1)), Aj'2 ^ ji (region A), a (0, 7r) columnar antiferro- 
magnet for 32 ^ (1 - 2ji)/(2(ji - 1)), Aj 2 > ji and 
ji < 1 (region C), a (7r, 0) columnar antiferromagnet for 
j 2 > (2ji - l)/(2(ji - 1)), Aji_> 1 (region B). In the 
remaining region D, the energy is minimized by two de- 
generate incommensurate spirals with pitch (Q x ,Q y ) = 
(2 arctan(a), 7r/2) and (27r — 2 arctan(a), 3ir/2), and a = 
Ui + V Ji + J2) /j2- Due to Qy's discretization all of the 
classical transitions are discontinuous. 

While the quantum analog of this rather rich phase 
diagram clearly deserves future analysis, our SE is con- 
fined to the region of < 1. Therefore we focus only 
this region regarding potential first order quantum phase 
transitions. First, the classical 'diagonal' transition from 
the (0,7r) to the (tt, it) state has no quantum analog, 
since it is connected to the region of <C 1. There, 
and instead of the classical (tt, 7t) or (0, tt) antiferromag- 
nets, the quantum model shows a phase of weakly cou- 
pled plaquette-singlets and the quantum ground state is 
protected by the singlet-triplet gap of the bare plaquette. 
To assess the relevance to the quantum case of the clas- 
sical transitions from the (0, tt) and (tt, tt) states into the 
doubly-degenerate (Q x ,ir/2), (Q x , 3tt/2) phase, we resort 
to a DMRG analysis of the ground state energy. For this 
we use the ALPS package^. Since we are only concerned 
with first order transitions, we refrain from any detailed 
finite size scaling analysis. We find that for ji^ < 2 there 
is negligible finite size dependence of the ground state en- 
ergy for FFSTs of lengths L w 30 ... 40, i.e. 120 .. . 160 
spins, and that m = 300 . . . 400 states kept in the density 
matrix lead to 4 ... 5 digits of precision, which is sufficient 
for our purpose. 

The inset in Fig. [2] shows a typical result for the ground 
state energy eo obtained from DMRG versus ji along the 
diagonal ji = j% ■ Obviously the ground state energy dis- 
plays a kink at ji = j° « 1, which we identify with a first 
order transition. The circles in Fig. [2] summarize a scan 
of locations of this transition which we have performed. 
The size of the circles is a rough measure of the numer- 
ical accuracy for the locations of the transition. These 
locations are remarkably close to those of the classical 
model. Below these transitions SE based on decoupled 
plaquettes is applicable. Beyond 312 ~ 1-5 the critical 
points are increasingly hard to detect accurately from the 
DMRG data. We speculate that below the first-order 
transition lines at if (,72) and J2O1) the bare plaquette 
state is adiabatically connected to a Luttinger liquid for 
32 — > 00 and ji — > 00, respectively. 



III. SERIES EXPANSION METHOD 

The main focus of this work is on SE in terms of ji and 
j 2 , starting from the limit of isolated plaquettes. To this 
end we decompose the Hamiltonian of the FFST into 

H = H + V, 
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Table I: Spectrum of the single plaquette. It consists of four 
equidistant energy levels E n — q„ — 2 labeled by the quantum 
number q n = 0, ..,3, the total spin S, and its z-component 
S z . The last column enumerates the states. 

H = £>o,«, V = V! + V 2 , (2) 

n 




Figure 3: Ground state energy per site eo versus j\ and j2, 
showing a monotonously decreasing behavior in the parameter 
range shown. Along the line of maximum frustration ji = , 
the energy gain is smallest. 



where /io,n is the plaquette Hamiltonian at site n 

ho, n = [Si • S 2 + S 2 • S 3 + S 3 • S 4 + S 4 • Si]„ (3) 

and the perturbation V = V\ + V2 is given by 

4 

Vi = ji Sj, n ■ Sj ;n +i, 

n 2—1 

^2 = j-2 / ](Sl |Tt • S2,n+1 + S2 : n • S3 ; „+i 

n 

+ 83,,! • S4 jn +i + S4 : „ • Si iTl +i 

)■ (4) 

The eigenstates of an isolated plaquette consist of four 
equidistant energy levels E n = q n — 2, labeled by the 
quantum number q n = 0, ..,3, and can be classified ac- 
cording to the total and the z-component of the plaque- 
tte spin S n = J2i=i Si,n- Table H] lists that the ground 
state is a singlet, the first excited state at q n = 1 is a 
triplet, the q n = 2 sector is composed of a singlet and 
two triplets, and for q n = 3 one quintet remains. 

In turn, Hq displays an equidistant spectrum, labeled 
by Q = Tl n qn- At V = 0, the Q = sector refers 
to the unperturbed singlet ground state |0) = ® n \so,n) 
composed of q n — singlets on all plaquettes. The 
Q = 1 sector comprises single |tg * m ) triplets inserted 
into |0) at site m. The Q > 2 sectors are of multi- 
particle nature. The perturbation V can be rewritten 
as V = Ylt=ijil2n=-N T n) where T„ represent raising 
(n > 0) or lowering (n < 0) operators within the spec- 
trum of H . For the FFST we find N = 4. 

It has been shown quite generally^, that models with 
the preceding type of spectrum allow for SE through a 
continuous unitary transformation (CUT) using the flow 
equation method of Wegner— . The basic idea is to trans- 
form H onto an effective Hamiltonian H e ff which is 
block- diagonal in the quantum number Q. This trans- 
formation can be achieved exactly order-by-order in j'12 
leading to 

H cfi = H Q + C»,mii -"*#*> (5) 

n,0<m<n 



where the C n<m are weighted products of the opera- 
tors which conserve the Q-number and have their weights 
determined by recursive differential equations, see~ for 
details. Due to Q-number conservation several observ- 
ables can be accessed directly from H e ff in terms of a 
SE in ji t 2- For systems with coupled spin-plaquette CUT 
SE has been used for one^S, two^Ir— and three^ dimen- 
sions. 



IV. RESULTS 

In this Section we present our findings from SE up 
to Q = 2, sectioning the discussion according to A. the 
ground-state energy, B. the one-, and C. the two-particle 
excitations. To assess the quality of the SE we com- 
plement our analysis by DMRG calculations for selected 
cases. 



A. Ground State Energy 

First we consider the ground state energy Eq. Q- 
conservation leads to 

E Q = (0|£T cff |0) , (6) 

where |0) is the unperturbed ground state. 

Evaluating this matrix element on chains with periodic 
boundary conditions (PBC) of a length L, sufficient not 
to allow for wrap-around of graphs with length N, i.e. 
L = N + 1, one can obtain analytic SEs which are valid 
to O(N) with respect to eqn. ([5]) in the thermodynamic 
limit. Evaluating the ground state energy per spin eo = 
E /(AL) up to 0(7) we get 

1 Ujf 17jf hh 25 Jr 25j| j\n 
C ° 2 96 96 3 384 384 16 
jig lglg# 1919j 2 4 _ 313^2 _ 313jij| 
16 96768 96768 1728 1728 
15457j^| 1510499j 1 5 1510499jf 39353j 4 j 2 
+ 48384 + 27095040 + 27095040 188160 
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39353j 1 j| , 4129273^^1 , 4129273j?j| 



188160 



27095040 



27095040 



522374480359 522374480359jf 92947333 jfja 



14748372172800 
92947333jijJ 



14748372172800 4267468800 
1494532466633j^| 



4267468800 4916124057600 
1494532466633ffi| 9888732599jf j\ 

4916124057600 
535161937582507j 1 7 



17069875200 
535161937582507jJ 



24777265250304000 24777265250304000 
3149204376698497jfj 2 3149204376698497jij| 



12388632625152000 
274272578154571 

412954420838400 
1186859862395537^^1 

2753029472256000 



12388632625152000 
274272578154571 j^f 
412954420838400 
1186859862395537.7^1 
h 2753029472256000 



(7) 



Here, the first term corresponds to the bare energy per 
spin listed in table HI Since the FFST can be mapped 
onto an identical FFST with j\ -H> 32 by a 7r/2-twist 
of the plaquettes around the tube, one expects that 
e o(iij 32) = eo(j2,Ji), which is obviously fulfilled. Figure 
[3] shows eo(0 < ji < 1,0 < 32 < 1) to be monotonously 
decreasing with ji^- Along the line j\ = j'2, the energy 
gain is smallest. Speaking differently, along its maxi- 
mally frustrated direction in parameter space the energy 
landscape is flattest. 

In Fig. [4] we asses the accuracy of eqn. in several 
ways. First, in panel (a) we compare eo with DMRG 
calculations along j 2 = where we expect no first-order 
transition. The DMRG results have been obtained for 
L = 30 and m = 300. Obviously the agreement between 
SE and DMRG is very good for ji < 0.75. In panel (b) we 
display a similar comparison along j\ = j2- Here the SE 
agrees very well with DMRG to even larger values of ji j2 , 
however only up to the first-order transition discussed 
in section [TTJ Third, in panel (a) we also include three 
DlogPade approximants to eqn. ([7]), which, similar to the 
DMRG, start to depart from the SE for j 1 > 0.75 and 
may be used to increase the interval of confidence of the 
SE only slightly by ~ 10% as evident from the DMRG 
results. For reference panel (a) also depicts the ground 
state energy per site of the ID HAFM with exchange 
constant j\ which is (— ln(2) + 1/4) ji. Evidently, in the 
parameter range we consider, the FFST is far away from 
the decoupled chain limit. 



B. One-particle Excitations 

Now we turn to the dispersion of one-particle eigen- 
states. By the latter we mean excited eigenstates of the 
effective Hamiltonian which, apart from having fixed Q, 
lattice momentum fc, and total spin quantum numbers 
5, to, are linear combinations of single entries of table 
U only. Q = 1 eigenstates are single particle states by 
construction. For Q 1, Q-conservation does not pro- 
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Figure 4: (a) Comparison between eo from plain SE (red solid 
lines) with DMRG for L = 30 and m = 300 (blue crosses), 
along ji = 0. Black dashed lines refer to three DlogPade 
approximants to eqn. 0. For reference the ground state en- 
ergy of the ID HAFM is also depicted (green dash-dotted), 
(b) Comparison of plain SE (red solid lines) with DMRG (blue 
crosses) along ji = ji. Up to the first-order transition dis- 
cussed in section [Til the agreement is excellent. 



tect single states in table U from decay into two states by 
virtue of H e ff. Eg. a \su) singlet with Q = 2 at site I 
could decay into two ^0^*)) triplets with total spin 5 = 
and Q = 1 + 1 at sites (see also appendix IVII[) . For 
the FFST however, and for Q = 2, S = 0, 1 and up to 
0(7), we find that all matrix elements of H e ff inducing 
such decay vanish identically. This feature can be traced 
back to the C4 symmetry of the tube. In fact, we obtain 
that changing e.g. the exchange couplings ji into j[ =/= j\ 
on one of the legs, renders the one-particle Q = 2 states 
unstable against decay. In summary, for each momentum 
k the Q — 2 spectrum contains three genuine one-particle 
levels with S = 0, 1, all of which are degenerate in to. For 
the remaining fwo-particle states with Q = 2 we refer to 
the next section. 

For the rest of this section we label the one-particle 
states by \i)^' Sn , where i refers to the plaquette's site 
and S n is the total spin, where the index n is only due 
to the fact that for Q = 2 there are two 5 = 1 states, 
say, n = a, b. Spin-z quantum numbers m are discarded 
because of SU(2) invariance. Due to Q-conservation and 
SU{2) invariance the effect of H c tf on \i)Q' Sn is limited 
to 



H off \3) C 



t 



Q,S n ,S„ 



\3+*) C 



(8) 



i,S n =S„ 



5 



which implies a shift in real space, and potentially a mix- 
ing of states of equal S with identical Q. The hopping 

OS s 

amplitudes t^' m do not depend on j due to trans- 
lation invariance. Therefore, by Fourier transformation 
|£)Q,S„ _ l/y/~LJ2j exp(-ikj)\j) Q - Sn we get the disper- 
sion from 



£^ pt o c (k) = Q ^(k\H eS \k} Q ' s " 



J Q,S„,S„ 



E S, 



s n ,s„ 



rQ,S n ,S„ 



t 







+ 2j2t? >3n ' Sm g(ik), (9) 



where, obviously tf' S "' Sm = t®f"' Sm for S n = S m . How- 
ever for S n 7^ S m , i.e. for the two Q = 2, S — 1 states, we 
find t?' S "' Sm = -t Q f"' Sm . In turn, g(x) is cos[i sin] (a;) 
for S n - [off] diagonal transitions. 

To obtain hopping amplitudes valid to O(N), in the 

OS s 

thermodynamic limit, the tj°' and the ground state 
energy E = (Q\H c g\Q) in eqn.flSJ) have to be evaluated 
on clusters with open boundary conditions (OBC), large 
enough to incorporate N-th order graphs for hopping pro- 
cesses of distance Depending on i, these clusters 
are of either of length N or N + 1. We have calculated 
analytic expressions^ for s (k) to 0(7) in ji^. 
For (Q,S) 7^ (2,1) eqn. © is already diagonal in S n , 
S m with eigenvalues E^g (k) = Eq^ s (k). Only for 
(Q, S) = (2, 1) eqn. © displays a 2 x 2-matrix struc- 
ture, referring to n = a, 6, with eigenvalues E]^\ a (k) and 

^Q P l b (fc)- 

In Fig. [5] we show the one-particle dispersions for se- 
lected values of j\ t 2- This figure displays only bare SE re- 
sults and no Pade extrapolations. Several comments are 
in order. First, the figure does not only display E^ s (k) 
for all S and Q < 2, but for curiosity also the quintet 
with Q = 3, assuming that the latter does not decay into 
multi-particle states - which we have not checked. 

As for the ground state energy, since exchanging j\ <-> 
j 2 maps the FFST onto an equivalent one by a 7r/2-twist 
there are related symmetries of the one-particle disper- 
sions. These respect the additional fact, that the single- 
particle states on the bare plaquette are the eigenstates 
of a four-site spin-1/2 chain with PBC, which carry one 
out of four momenta fc n j_ = nir/2 with n = 0,1,2,3 
transverse to the FFST. In turn, exchanging j\ <R- j2 
maps the one-particle dispersions onto identical ones up 
to a shift of the Brillouin zone by one out of k n ±, and 
may also exchange the dispersion branches for the de- 
generate bare Q, S = 2,1 states. We have checked 
this to be fulfilled by all E^* s Jk). Eg., for (Q,S) = 

(1,1) we have E^* Sn (k,ji,j 2 ) = E^Jk + n,j 2 ,ji) 
and for (Q,S) = (2,la[6]) we have E^ la[b] (k,j 1 ,j 2 ) = 
E QMa] (fc+H 7r /2,J 2 ,ii ) . Additionally E^ la[b] (k,h, j 2 ) 
are degenerate at j 2 = 0. 

Considering the elementary triplet dispersion E\ p ^(k) 
at j 2 = 0, we observe a zone-boundary gap in Fig. [5] (a) 
which decreases as j\ increases, see panel (b). For four- 
leg ladders without frustration, i.e. j 2 = 0, the analytic 
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Figure 5: Genuine one-particle dispersions on the FFST for 
Q < 2 at various values of (j'1,,7'2). These comprise three 
triplets (red circle, green triangle up, green triangle down) 
and one singlet (blue diamond). The Q=3-quintet (magenta 
x) has not been tested against multi-particle decay and is 
only shown for reference. DMRG results at zone-boundaries 
are shown with large crosses. 



expression for the gap we get is 



P ipt/ ■ nA , 4 41 . 2 349 , 3 
4596401 . 160407997 



4596401 . 4 169497997 , 5 

^39191040 : ' 1 ~ 4702924800 Jl 
689874137639377 . 6 

~2986545364992000 Jl 
8430165345498432721 , 7 

~ 45156565918679040000 Jl ' 



(10) 



which, in passing, improves earlier SE to 0(4) on un- 
frustrated four-leg ladderJ*^ by three orders. The disper- 
sion E\ p ^{k) shows rather strong oscillations in panel (b). 
These are not related to convergence issues of the SE, but 
are robust features of hopping to more than only nearest- 
neighbors. 

Increasing j 2 from zero, at fixed j\, as from Fig. 0(a) 
to (c) the bandwidth of El P *(k) is reduced. On the line 
of maximum frustration, i.e. at ji = j 2 in panel (d), it is 
only very small, albeit not exactly zero, namely O(jf). 
A similar tendency can be observed in all other (Q, S) 
sectors, yet less pronounced. 

To assess the quality of the SE, we have also calculated 
the elementary (Q, S) = (1, 1) gap by DMRG. The results 
are shown by the large crosses at the zone-boundaries in 
Fig. [S] As is obvious, the agreement is very good in all 
four panels. This is particularly noteworthy for panel (b), 
corroborating the preceding statement, that the oscilla- 
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Table II: Q = 2 states from the bare plaquette. Cf r ^ l lm _, I = 
(In, lm — nlSm) are Clebsch-Gordan coefficients. 

tions in ^^(/c) for larger values of j2 are robust features 
of the SE and no convergence issue. 

C. Q=2 Two Particle States 

In this section we focus on two-particle states gener- 
ated from two one- particle triplets out of the (Q,S) = 
(1,1) sector. Apart from generating a continuum of 
states at fixed total momentum with respect to the rela- 
tive momentum, interactions may lead to additional col- 
lective (anti) bound states which split off from the contin- 
uum. The existence of such collective states is know for 
simpler ID quantum magnets, such as chains and two-leg 
ladders^— . Here we show that similar states exist also 
on the FFST. To this end we evaluate the two-particle 
spectrum following ideas developed for dimer-SE of two- 
leg ladders 3 -^. For the sake of completeness technical de- 
tails of this procedure are revisited in appendix IVHI 

On the FFST, the Q = 2 sector allows for one-particle 
and two-particle states. This is different from dimer sys- 
tems where only two-particle states occur for Q = 3^. 
The real-space representations of all Q = 2 states with 
proper spin quantum numbers S, to are listed in table 
HIl Here the genuine two-particle states are labeled by 
\hj) Sm > where i,j refers to sites on the ID lattice and 
S, to to total spin, and spin-z quantum numbers. As dis- 
cussed in the preceding section the one- and two-particle 
excitations at Q = 2 do not mix. Therefore we focus on 
the genuine two-particle excitations hereafter. 

Because of translational invariance the two-particle 
states can be classified according to a center-of-mass mo- 
mentum q and a relative distance d 

\q,d) Sm = ^=J2 eiq{r+d/2) \ r ' r + d ) Sm ( U ) 

These states have an exchange parity P = (— l) s for 
d -> -d, i.e. \q,-d) Sm = P\q,d) Sm . The effective 
Hamiltonian, eqn. ([5]), will mix states at different d. The 
two-particle spectrum is obtained from the sum of the 
one- and two-particle irreducible Hamiltonians (se o 30 i 32 
and appendix IVII|) . labeled by Hi and H2, with analytic 
matrix elements calculated by SE from the states of eqn. 
([IT]) . We have evaluated these matrix elements to 0(6). 
For each q, this leads to an eigenvalue problem with a 
rank set by the number of lattice sites (or relative mo- 
menta) L. This can be diagonalized numerically on finite, 
but large lattices. See appendix IVHI for definitions and 
formal details. 
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Figure 6: Two-particle density of states close to the zone 
boundary at q = 99tt/I00 on an L = 200 site FFST for 
(jijjz) — (0.6,0.4). Delta functions at each two-particle en- 
ergy broadened by 77. (a) reducible spectrum for S=0 and 
S=l . (b) difference between bare and reducible spectrum for 
S=0 and S=I. (c) complete spectra for S=0,f ,2. 



In Fig. [6] we show the two-particle spectrum ps(q, <jS) = 
^2 p S(u! — Eg^fapy) at fixed total momentum q and 
spin S, to, summed over the relative momentum p, where 
Eg P ^(q,p) are the two-particle eigenenergies. We choose 
(J11J2) — (0.6,0.4) and q ~ n as an example. To ana- 
lyze the effects of the interaction we consider Ej'^ n (q,p) 
at various levels of approximation. First, we set the ir- 
reducible two-particle interactions H2 to zero - i.e. no 
'actual' interactions occur. The corresponding spectrum 
is shown in Fig. HI (a) for S — and 5 = 1. These spec- 
tra are indistinguishable on the scale of this plot. The 
van-Hove singularities are related to the extrema of the 
one-particle dispersion in the (Q, S) = (1, 1) sector as in 
Fig. [5] Even though the irreducible one-particle interac- 
tion Hi will act only on one of the two particles, this does 
not imply, that the spectra in Fig. [S] (a) are identical to 
the convoluted bare spectra of two single |£ *)-particles 
at fixed total momentum q and finite L. This is due 
to the hard-core constraint which forbids double occu- 
pation of plaquettes by |iQ~)-particles. This constraint is 
encoded in the matrix elements of Hi when acting on the 
two particle sector. Since there are L double occupancies 
within L 2 two-particle states, their removal amounts to 
a 1/1 effect on the integrated spectral weight at fixed to- 
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Figure 7: Diagonal element of the irreducible interaction 
{q, d\H2\q, d) versus relative distance of two particles for S = 
0, 1, and 2, and for q = and ir, multiplied by 10 d_1 . (The 
matrix elements of Hi at 0(6) are zero for d = 5, 6.) 



tal momentum q. At fixed total and relative momentum, 
this effect will be quantitatively identical for two-particle 
states of identical exchange parity (— ) s . For this reason 
Figs. [6] (a) and (b) contain only S = and 1 sates. Fig- 
ure |6] (b) depicts the difference between the bare spectra 
and that obtained from Hi for £ = and S = 1 for 
L = 200. Clearly the integrated spectral weight for both 
values of S is finite but of 0(1/ L) at most. Differences 
are visible between the S = and S = 1 spectrum - al- 
beit very small. Figure [6] (a) shows that no (anti-)bound 
state arises in the two-particle sector solely due to the 
hard-core repulsion. 

Next we consider the complete Q = 2 spectrum in Fig. 
[6] (c) for S = 1,2 and 3, including the irreducible two- 
particle interaction H2. The latter is spin- dependent and 
longer-ranged. This figure clearly demonstrates a main 
point of this section, i.e. the occurrence of three split- 
off (anti)-bound states (above) below the continuum de- 
pending on their spin. Apart from that, and similar to 
the hard-core constraint H2 leads to an 0(l/L) redis- 
tribution within the two-particle continuum. The total 
momentum q in Fig. [5] has been chosen, such as to ev- 
idence the (anti)-bound states clearly. Depending on q 
they can approach and also merge with the continuum. 

A qualitative understanding of binding versus anti- 
binding of the collective two-particle states can be ob- 
tained from considering the diagonal matrix elements 
Sm (q, d\ H2\q, d) s m of two particles versus their relative 
distance d at fixed total momentum q. This matrix ele- 
ment is oc Ss's^m'm and independent of m. It is shown in 
Fig. [7]for ji t 2 identical to that of Fig. [6l First, this figure 
demonstrates that the irreducible two-body interactions 
are very rapidly decaying as a function of the relative 
distance of the particles (note the scaling of the matrix 
elements by 10 d_1 on the y-axis). Second, and for both 
total momenta depicted, the interactions for S — 1 are 
attractive, while those for 5 = 0, 2 are repulsive. This is 
consistent with binding versus anti-binding of the collec- 
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Figure 8: Momentum resolved spectrum of all states on the 
FFST with Q < 2 for (ji,j 2 ) = (0.6,0.4) and (0.6,0). Apart 
from four genuine one-particle states, i.e. three triplets (black 
dashed, blue solid, green solid), and one singlet (black solid), 
there are three two-particle continua with S=0,l,2, and three 
collective two-particle states: a bound triplet (magenta dots), 
an anti-bound singlet (black dots) and an anti-bound quin- 
tet (red dots). Only ~ f0% of the scattering states between 
the lower and upper bounds of the continua are depicted to 
avoid graphic cluttering. DMRG results for one and two- 
particle states at zone-boundaries(centers) are depicted with 
large crosses. 



tive states in Fig. [5] Yet, one should keep in mind the 
matrix elements s m (q 1 d! 7^ d\ i?2|s, d) Sm which are off- 
diagonal in d' ,d. We have not tested if these invalidate 
the simplified argument given here. Finally, we empha- 
size, that the ordering of anti-bound versus bound states 
on the FFST is different from that on plain two-leg lad- 
ders. In the latter case two-particle singlets also form 
bound states with a binding energy larger than that of 
the two-particle triplets. 

In Fig. [8] we show the dispersion of all one and two- 
particle states from the Q = 1 and Q = 2 sector as 
a function of total momentum for two sets of j\ t 2, ob- 
tained by SE, together with DMRG results for one and 
two-particle states at the zone-boundary (center). This 
figure is another main result of our paper and summa- 
rizes several aspects. First, the spectrum is very rich 
and consists of various discrete states and three superim- 
posed two-particle continua. This should be contrasted 
against plain two-leg ladders, which show a less involved 
low-energy spectru m 28 i 32 . Some of the discrete states are 
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genuine one-particle states and some are collective two- 
particle states. For the parameters depicted, the collec- 
tive two-particle states are clearly visible only close to 
the zone boundary. The potential existence of critical 
parameters jf 2 or wave- vectors q c for which the collec- 
tive states merge with the continua is unclear at present. 
The Q = 2 one-particle states are contained almost com- 
pletely within the spectral range of the continua and 
split off from the latter only in a range of q similar to 
the collective states. This situation is also very different 
from two-leg spin-ladders in which the states which split 
off from the continuum close to the zone boundary are 
collective (anti)-bound states only. From an experimen- 
tal point of view this may pose a challenge on discrimi- 
nating between such genuine one-particle and collective 
two-particle states. Finally, the effects of frustration are 
clearly visible in going from Fig. [8] a) to b): the com- 
plete spectrum shows a tendency to localize. However, in 
contrast to the one-particle dispersions the two-particle 
interactions are not reduced approaching the line of max- 
imum frustration, ji = ji- In turn the relative splitting 
of the collective states from the continuum is enhanced 
by frustration. 

To conclude this section we provide some measure of 
convergence of the SE for the two-particle continuum in 
Fig. [U This figure shows the relative change in energy of 
all two-particle states at a fixed total momentum versus 
their energy when switching from an 0(6) to an 0(7) 
evaluation of the matrix elements of H\ 2- For this we 
confine ourselves to S = 0. As is obvious, the changes 
are completely negligible and fully justify the use of 0(6) 
SE for the two-particle states. We note in passing, that 
0(7) SE will not only lead to longer-range irreducible 
one-particle hoppings, but also increase the range of the 
irreducible two-particle interactions in Fig. retaining 
however the rapid decay with distance depicted there. 
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Figure 9: Relative error in all two-particle energies for an 
FFST of 200 plaquettes at fixed total momentum q for switch- 
ing from 0(6) to 0(7) SE. (The apparent 'two-curve' struc- 
ture is due to an oscillation of the error.) 



sidering the weak rung-coupling limit. 
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V. CONCLUSIONS 

To summarize, we have studied the strong rung- 
coupling regime of the frustrated four-spin tube. Re- 
markably, this spin model displays several twists with 
respect to conventional two-leg spin ladders. It allows 
for frustration induced quantum-phase transitions into 
phases other than the plaquette singlet phase. The na- 
ture of these phases remains to be explored. The struc- 
ture of the excitations is both richer, and different as com- 
pared to two-leg ladders, including a reordering of two- 
particle (anti)bound states and additional elementary ex- 
citation modes, some of which can be made to decay by 
reducing the symmetry of the FFST. Finally geometric 
frustration on the FFST introduces a control parameter 
allowing for almost localization of excitations and a flat- 
tening of the energy landscape. As compared to potential 
material realizations of the FFST in Cu2Cl4-DgC4S02, 
several extensions have to be considered in the future. 
These comprise e.g. reducing the C4-symmetry and con- 



VII. EFFECTIVE HAMILTONIAN MATRIX IN 
THE Q = 2 SECTOR 

In this appendix we list some of the details necessary to 
evaluate the two-particle matrix elements by SE. First, 
we refer to the parity of the states in eqn. (fTTj) . Be- 
cause of this we may confine ourselves to d > (d = 0, 
i.e. \r,r) is forbidden). The matrix elements of 
the one(two)-particle irreducible Hamiltonian i?i(2) are 
obtained by subtracting the one(two)-particle reducible 
contributions from the matrix elements of H e f f ' 

to,n = *n,o = (n|JTi|0) = (n\H eff - H o \0) (12) 

a 0.n ~ ^0,nEf 

t d \n,d> = (n,n + d'\H 2 \0,d) (13) 
= (n,n + d'\H eff - H x - H Q \0,d) 

= a d\n,d> ~ 8o,n8d,d' E ~ $d,n+d'to n - 
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(n\H eff \0), a* nd , = (n,n+d'\H eff \0,d) and 
cl refers to the cluster on which these matrix elements are 
evaluated. For io,r» and td\ n ,d' to be size consistent, the 
cluster has to be the largest linked cluster for a given one- 
and two-particle state in real space at a given order N of 
the SE. The ig' n in eqn. (fT3")l refers to t .n evaluated on 
the same linked cluster 'cl' as the a^j n d , corresponding 



to a given td\ n d' and Pq 



\Hq\0) is the ground state 



energy of that cluster. The two-particle spectrum results 
from the eigensystem of 



(q,d'\H eff -H \q,d) 1 



(q,d'\H2 + Hi\q,d) s (14) 



which, due to translational and spin-rotational invariance 
is diagonal in q,S, and m and is independent of m. For 
each (q, 5*, m) the two matrices on the left of eqn. (fH)) 
are hermitian with, matrix-indices (d',d). Thus, for the 
remainder we may consider only d' ^ d. We begin with 
hid'd = S (q, d'\H2\q, d) s . Because a linked cluster can 
have at most N + 1 sites, fi2d'd = for d',d > N. I.e. 
^2d'd is an N x N matrix. The action of H2 is 



H 2 \q,d) S = -Lj2 



e iq{r+d/2) x 



E td\n,d'\r + n,r + n + d') s 

E t dln , d ,e^ d - d '^\q,d') s (15) 

max(7i+d',rf- n)^-N 

here n must be restricted to max(n + d',d— n) ^ iV, 
because of the sites (0,d, n, n + <i') to reside within a 
linked cluster at order N . I.e. n is confined to the interval 
d — N ^ n ^ N — d' , which has its midpoint at n m = 
(d-d')/2 



H 2 \q,d) S = { 

ir 

n=(d-d')/2e 

2/ 

iV 

n>(d-d')/2eZ 



E td\(d-d')/2.d'- 
max(n+d ,d—n)^.N 



(16) 



max(n+d',d-n)^iV 



cos[<7(^-^-n)]}M') S 



where td\ n ,d' = ^Id-d'-n.d' has been used, which refers to 
reflection symmetry of the two-particle matrix element 
about the midpoint, and the addends with matrix ele- 
ments id|((2-d')/2,d' occurs only if n m = (d — d')/2 £ Z. 

Next, we apply Hi to a two-particle state. This will 
only shift one of the particles, i.e. 



Hi\q,d) S = J2 t£n( e ~ 19n/2 + el " n/2 ) 
n^d 

x _jL Y e Mr+(d-n)/2) t rr + d _ n \S 
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Figure 10: Linked clusters with size l c i and with initial (half 
circles) and final (wedges) two-particle states for all addend 
in eqn. (THJ at order N = 4 and for d = d' = 1. (I)R labels 
(ir)reducible graphs. 



E 

n^d 



,nq 



K.n cos(-^)sgn(d - |?, |d - n\) L 



(17) 



For each q, Hi connects states d and \d — n\, V — N ^ 
?i $J N. This is equivalent to a band-matrix of width 
2iV + 1. For d ^ N, the non-zero content of the columns 
of this matrix is independent of d, for 1 ^ d < N this 
is not so. The latter relates to the exclusion of on-site 
double-occupation. 

For reference we explicitly display one specific low- 
order matrix element s (q,d'\Hi + H2\q,d) s from eqns. 
PI OH EJ), namely d' = d = 1, N = 4 and S = 0. 
Figure [10] relates to the two-particle irreducible matrix 
element from H2, and depicts all td| n ,d' on their respec- 
tive clusters including a label for reducing contributions 
to be evaluated according to eqns. (fl"2| [T3|) . These are 
E§ , 4q 0, tf \ for n = and £§' 2 for n = 1 for this particu- 
lar matrix element. Note that on finite clusters tjq'q and 
t\ j are not necessarily identical. The one-particle irre- 
ducible matrix element from Hi is straightforward. The 
complete matrix element reads 



°(q,l\Hi+H 2 \q,l) =2 
18511597$ 



431$ 



5827$ 



156764160 



J2 

6 



+ : 1± 

6 ^ 432 ^ 15552 
461j lj2 5725$j 2 



216 



15552 



6030787$j 2 , 431$ 5725j lj2 2 



5598720 432 
51190589j 2 j| 5827j| 



26127360 
18511597$ 
156764160 



15552 
14jf 
27 



15552 

6030787j lj2 3 

5598720 
68j? , 72187J! 4 



243 



349920 



ggiijg 68j 2 j 2 23057j|i2 14j| 68ju 2 2 

27 243 17496 27 243 

21559j 2 i| 68jf 23057j lj2 3 , 72187$ 



9720 243 17496 349920 
{76jf 493j 3 j 2 , 278i 2 i 2 493jij| 



cos(q) 



V 729 



729 



243 



729 
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cos{2q) (18) 

Finally, in the Q = 2 sector the single-particle states 
l*2?)> an( ^ l s i«) nee d to be considered. Any one- 
particle irreducible Hamiltonian matrix element between 
the two- and one-particle states for Q = 2 is zero due 
to orthogonality. Two-particle irreducible contributions 
require matrix elements of the form 

Sm (q,d'\H eff \q) Sm , (19) 



where \q) Sm = ^ r e lqr \r) Sm / \/L is a one-particle state 
of spin quantum numbers S, m with Q — 2. This implies 
a decay of Q — 2 single-particle states into Q = 2 two- 
particle states. Up to the order that we have performed 
the SE, we have not observed such decay. In turn the 
Hamiltonian matrix in the Q = 2 sector of the four-tube 
is already diagonal w.r.t. the particle number. 
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